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Abstract 


Many  heavy  metals,  including  nickel  (Ni),  cadmium  (Cd),  and  chromium  (Cr)  are  toxic 
industrial  chemicals  with  an  exposure  risk  in  both  occupational  and  environmental  settings  that 
may  cause  harmful  outcomes.  While  these  substances  are  known  to  produce  adverse  health 
effects  leading  to  disease  or  health  problems,  the  detailed  mechanisms  remain  unclear.  To 
elucidate  the  processes  involved  in  the  of  toxicity  of  nickel,  cadmium,  and  chromium  at  the 
molecular  level  and  to  perform  a  comparative  analysis,  H4-II-E-C3  rat  liver-derived  cell  lines 
were  treated  with  soluble  salts  of  each  metal  using  concentrations  derived  from  viability  assays 
and  gene  expression  patterns  were  determined  with  DNA  microarrays. 

We  identified  both  common  and  unique  biological  responses  to  exposure  to  the  three 
metals.  Nickel,  cadmium,  chromium  all  induced  oxidative  stress  with  similar  genes  and 
pathways  responding  to  this  stress.  Although  all  three  metals  are  known  to  be  geno toxic, 
evidence  for  DNA  damage  in  our  study  only  exists  in  response  to  chromium.  Nickel  induced  a 
hypoxic  response  as  well  as  inducing  genes  involved  in  chromatin  structure,  perhaps  by 
replacing  iron  in  key  proteins.  Cadmium  distinctly  perturbed  genes  related  to  endoplasmic 
reticulum  stress  and  invoked  the  unfolded  protein  response  leading  to  apoptosis.  With  these 
studies,  we  have  completed  the  first  gene  expression  comparative  analysis  of  nickel,  cadmium, 
and  chromium  in  H4-II-E-C3  cells. 
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Introduction 


Many  heavy  metals,  including  nickel,  chromium,  and  cadmium,  are  widely  distributed, 
posing  occupational  and  environmental  exposure  risks  which  may  result  in  adverse  health 
effects.  Exposure  to  these  metals  can  occur  through  contact  with  contaminated  soil,  air,  water, 
and  food,  or  by  absorption  through  the  skin  as  a  result  of  manufacturing,  pharmaceutical,  or 
industrial  processes  or  environmental  contamination.  While  these  metals  are  known  to  cause 
adverse  health  effects  and  to  be  toxic  to  the  lungs,  kidneys,  liver,  and  other  vital  organs  [1,2],  the 
genes  and  toxicity  pathways  that  respond  to  metal  exposure  are  not  completely  known.  To 
further  elucidate  common  and  unique  mechanisms  of  toxicity  and  identify  the  genes  involved  in 
the  perturbed  pathways  we  performed  head-to-head  comparisons  of  the  gene  expression  profiles 
and  toxic  mechanisms  of  the  rat  liver  derived  cell  line  H4-II-E-C3  exposed  to  nickel,  cadmium, 
or  chromium.  We  identified  both  unique  and  common  biological  responses  to  metal  exposure  in 
this  comparative  analysis,  including  a  common  induction  of  oxidative  stress. 

Nickel,  cadmium,  and  chromium  have  been  used  in  a  wide  variety  of  products  and 
industrial  processes,  resulting  in  a  high  likelihood  of  occupational  and  environmental  exposure. 
Nickel  is  used  extensively  in  many  industrial  and  consumer  products  such  as  stainless  steel, 
magnets,  coins,  and  alloys;  evidenced  by  the  fact  that  882  of  the  1,662  current  sites  on  the 
National  Priorities  List  targeted  for  federal  clean-up  activities  contain  nickel  [3],  Chromium  is 
extensively  used  for  stainless  steel  production,  chrome  plating,  and  pigments  and  is  responsible 
for  500,000  industrial  exposures  in  the  United  States  [4,5],  Exposure  to  cadmium  can  occur  as  a 
result  of  mining,  metal  processing,  welding,  burning  fuels,  the  production  and  use  of  phosphate 
fertilizers,  and  leaching  of  metal  waste,  yet  tobacco  smoke  and  food  are  still  the  main  sources  of 
intake  [6]. 
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While  many  of  the  adverse  health  effects  of  nickel,  cadmium,  and  chromium  are  similar, 
the  exact  mechanisms,  modes  of  action,  and  biochemical  pathways  affected  by  each  metal  differ. 
For  example,  all  three  metals  induce  oxidative  stress,  but  nickel  and  chromium  undergo  Fenton 
type  reactions  fonning  reactive  oxygen  species  while  cadmium  is  thought  to  cause  oxidative 
stress  through  the  inhibition  of  antioxidant  enzymes  [7,8].  Similarly,  all  three  metals  have  been 
shown  to  be  genotoxic,  but  chromium  is  the  only  one  of  the  three  metals  shown  to  interact 
directly  with  DNA,  forming  Cr-DNA  adducts  and  causing  DNA  damage.  Nickel  and  cadmium 
are  thought  to  damage  DNA  through  the  inhibition  of  repair  enzymes  [7].  Nickel  and  cadmium 
deregulate  cell  proliferation  by  perturbing  various  signaling  pathways  and  transcription  factors, 
possibly  through  reactive  oxygen  species,  although  the  activation  of  these  pathways  is  poorly 
understood  [7]. 

In  order  to  identify  the  unique  and  common  effects  of  nickel,  cadmium,  and  chromium  at 
the  molecular  level,  we  perfonned  side-by-side  comparisons  of  the  effects  of  the  metals  in  H4-II- 
E-C3  cells  using  Affymetrix  DNA  microarrays.  H4-II-E-C3  cells  were  selected  for  use  as  they 
are  well  characterized  and  metabolically  active  liver  models  [9].  The  cells  were  exposed  to  nickel 
(II)  chloride  (NiCE),  cadmium  chloride  (CdCE),  or  sodium  dichromate  (Na2Cr207).  We 
identified  992  probe  sets  whose  expression  is  affected  by  exposure  to  at  least  one  of  the  metals 
(430  in  nickel,  456  in  chromium,  and  288  in  cadmium).  In  the  comparison  study  here  we 
demonstrated  that  the  metals  were  able  to  elicit  distinct  changes  in  the  gene  expression  profiles, 
and  we  identified  both  common  and  unique  mechanisms  of  toxicity  among  the  metals. 
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MATERIALS  AND  METHODS 


Cell  Culture  Conditions  and  Exposures 

H4-II-E-C3  cells  (ATCC,  Manassas,  VA)  were  grown  in  Dulbecco's  Modified  Eagle’s 
Medium  (DMEM)  (Lonza,Walkersville,  MD)  containing  10%  fetal  bovine  serum  (Invitrogen, 
Carlsbad,  CA)  and  10  mL  Glutamax  (Invitrogen)  in  T75  flasks  incubated  at  37  °C  with  5% 
carbon  dioxide.  Exposures  were  initiated  once  flasks  were  90  ±  10%  confluent  using  the  test 
chemicals  NiCfi,  CdCE,  and  Na2CriC>7  (Sigma-Aldrich,  St.  Louis,  MO).  Exposure  concentrations 
were  chosen  based  on  the  CellTiter-Fluor  Cell  Viability  and  CellTiter  96  Aqueous  Non- 
Radioactive  Cell  Proliferation  Assays  (Promega,  Madison,  WI)  at  the  no  observed  cell  death 
level,  and  at  the  20%  and  50%  cell  death  levels  corresponding  to  40,  140,  and  400  pM  for  NiCE; 
0.275,  1,  and  10  pM  for  Na2Cr2C>7;  and  0.2,  0.55,  and  1.2  pM  for  CdCE.  Prior  to  exposure  flasks 
were  washed  twice  with  serum  free  DMEM  to  remove  residual  serum  components  with  a  five 
minute  incubation  between  washes.  Fifteen  milliliters  of  serum  free  DMEM  containing  the 
proper  concentration  of  toxicant  were  then  added  to  each  flask  for  24  hours.  Four  biological 
replicates  were  performed  for  each  condition,  including  an  unexposed  control. 

RNA  Extraction 

The  cells  were  scraped  from  the  surface  of  the  flasks  and  were  homogenized  using  a 
Dounce  homogenizer  in  Trizol  solution  (Invitrogen).  Total  RNA  was  extracted  using  Trizol 
solution  per  the  manufacturer’s  instructions.  An  RNeasy  Midi  Kit  cleanup  (Qiagen, 

Gennantown,  MD)  was  performed  per  the  manufacturer’s  instructions  to  remove  residual  salts 
and  organic  solvents.  RNA  quality  and  quantity  were  detennined  using  the  Agilent  Bioanalyzer 
Series  II  RNA  6000  Nano  LabChip  Kit  and  2100  Bioanalyzer  (Agilent,  Palo  Alto,  CA). 
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Microarray  Preparation  and  Processing 


cDNA  and  labeled  cRNA  were  prepared  using  the  Affymetrix  GeneChip®  Two-Cycle 
Target  Labeling  kit  and  7.5  pg  total  RNA  according  to  the  GeneChip  Expression  Analysis 
Technical  Manual  (701021  Rev.  5).  Twenty  micrograms  of  biotin-labeled  cRNA  was  sent  to  the 
laboratory  of  Dr.  Maryanne  Vahey  at  the  Walter  Reed  Army  Institute  of  Research  Vaccine 
Genomics  Laboratory  for  processing  and  scanning  on  the  GeneChip  Rat  Genome  230  2.0  Array 
using  Affymetrix  instrumentation  according  to  the  manufacturer’s  instructions  (Affymetrix, 
Santa  Clara,  CA). 

Data  Analysis 

Microarray  data  was  processed  for  background  adjustment,  normalization,  and 
summarization  using  the  Robust  Multi -Array  Averaging  method  (RMA)  [10]  using  Partek 
Genomic  Suite  software  (Version  6.4  Copyright  2009,  St.  Louis,  MO).  Raw  data  files  can  be 
found  in  the  NCBI  Gene  Expression  Omnibus  (accession  number  to  be  detennined).  The 
microarray  data  was  examined  for  outliers  using  a  principal  component  analysis  (PCA)  in  Partek 
6.  Pairwise  correlation  analysis  and  inter-replicate  dot  plots  of  all  probe  sets  were  performed  to 
verify  reproducibility.  Replicates  were  accepted  with  an  R  >0.95  and  no  gross  deviations  from 
linearity.  If  a  sample  did  not  meet  these  criteria,  a  new  microarray  was  processed  from  the  total 
RNA.  A  present,  absent,  or  marginal  detection  call  for  each  probe  set  was  determined  using  the 
Affymetrix  GCOS  algorithm,  and  only  probe  sets  with  a  present  detection  call  for  all  samples  in 
at  least  one  condition  were  retained  for  analysis  [11]. 

An  analysis  of  variance  (ANOVA)  was  perfonned  to  detennine  which  genes  were 
differentially  expressed  due  to  treatment.  The  16,026  probe  sets  that  met  the  present  detection 
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call  criteria  were  analyzed  using  2-way  ANOVAs  (dose  and  batch)  with  contrasts  for  each 
exposure  concentration  versus  the  control  using  Partek  for  each  metal.  The  batch  variable  was 
included  to  control  for  differences  observed  in  the  PCA  resulting  from  different  experimental 
and  processing  dates.  Probe  sets  with  a  Benjamini  and  Hochberg  False  Discovery  Rate  (FDR) 
[12]  less  than  or  equal  to  0.001  for  the  concentration  variable  and  a  1.8  or  greater  fold  change 
from  control  in  at  least  one  treatment  condition  were  retained  for  bioinformatic  analysis. 

After  an  initial  unsupervised  ontology  analysis,  it  was  observed  that  multiple  related 
categories  were  present  in  the  results  that  were  similar  to  known  effects  of  these  metals. 
Therefore,  a  manual  binning  method  was  devised  in  order  to  attribute  intuitive  biological 
functions  to  a  large  portion  of  differentially  expressed  genes.  This  scheme  assigned  the  major 
biological  processes  that  were  modulated  by  treatment  with  the  toxicants  by  sorting  probe  sets 
into  groups,  or  bins,  based  on  Gene  Ontology  (GO)  categories  and  the  known  effects  of  the 
metals.  Seven  bin  categories  were  created:  cell  cycle,  oxidative  stress,  ion  homeostasis, 
apoptosis,  energy  regulation,  hypoxic  response,  and  DNA  damage,  replication,  and  repair.  Each 
bin  was  compromised  of  multiple,  related  GO  terms  based  on  the  GO  biological  process  terms 
provided  by  Affymetrix  in  the  annotation  file  (build  29,  2009-7-13)  for  the  Rat  Genome  230  2.0 
Array.  The  GO  terms  found  in  each  bin  can  be  found  in  the  supplemental  data  (Table  SI).  Probe 
sets  were  assigned  to  a  bin  if  the  GO  tenn  associated  with  that  probe  set  was  also  contained  in 
that  particular  bin.  A  chi-squared  test  was  used  to  test  bin  enrichment  (p< 0.05).  Probe  sets  that 
did  not  contain  any  biological  process  annotation  were  not  considered  for  significance  testing. 

Differentially  expressed  probe  sets  were  clustered  using  Vxlnsight  and  VxArrayhnport 
0.2.5  with  default  settings  [13]  (Sandia  National  Laboratories,  Albuquerque,  NM)  to  identify 
probe  sets  with  similar  expression  patterns.  Vxlnsight  uses  a  force  directed  placement  algorithm 
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to  move  similar  items  closer  together  while  simultaneously  pushing  dissimilar  objects  away  from 
each  other,  and  then  displays  the  relationships  on  a  3D  terrain-like  map  [14].  Clusters  were 
manually  selected  by  their  natural  boundaries  using  the  terrain  view.  Transcription  factor 
enrichment  for  probe  sets  in  each  cluster  and  metal  were  investigated  using  Metacore  (GeneGo, 
St.  Joseph,  MI). 

Ingenuity  Pathway  Analysis  (IP A)  software  (Ingenuity  Systems,  www.ingenuity.com, 
analysis  date  2009-1 1-09)  was  used  to  explore  the  biological  implications  of  the  data.  Core 
analysis  were  perfonned  on  the  data  using  the  Rat  Genome  230  2.0  Array  as  the  reference  set 
with  all  other  default  settings  selected.  We  considered  networks,  canonical  pathways,  and 
biological  functions  statistically  significant  with  a  />value<0.05.  A  canonical  pathway  in  IPA  is  a 
well  characterized  metabolic  or  cell  signaling  pathway  that  was  drawn  based  on  the  IPA 
Knowledge  Base  while  bio  functions  associate  biological  functions  and  complex  interactions 
from  the  Ingenuity  Knowledge  Base  with  the  experimental  data  set. 

RESULTS  AND  DISCUSSION 

Since  nickel,  cadmium,  and  chromium  are  potential  environmental  and  occupational 
hazards,  we  undertook  a  study  to  identify  common  and  unique  mechanisms  of  toxicity  for  the 
three  metals  with  a  focus  at  the  level  of  gene  expression  and  molecular  pathways.  We  exposed  a 
rat  hepatoma  derived  cell  line  (H4-II-E-C3)  to  three  concentrations  of  NiCb,  CdCh,  or  Na2CriC>7 
for  24  h  and  then  analyzed  for  transcriptional  changes  using  whole  genome  DNA  oligonucleotide 
microarrays.  Pathways  and  biological  functions  affected  by  the  exposure  to  each  metal  were 
indentified  and  then  compared  among  the  metals  to  further  explore  similarities  and  differences  in 
the  responses  to  the  three  metals. 
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Microarray  Analysis 


To  identify  genes  differentially  expressed  due  to  exposure  to  the  metals,  we  measured 
mRNA  levels  using  whole-genome,  DNA  oligonucleotide  microarrays.  The  data  was 
preprocessed  using  the  RMA  method  and  filtered  to  select  only  probe  sets  with  a  present  call  in 
all  replicates  of  at  least  one  condition;  16,3 1 1  out  of  the  3 1,099  possible  probe  sets  were  retained 
for  further  analysis.  Differentially  expressed  genes  were  identified  by  calculating  two-way 
ANOVAs  (dose  and  batch)  for  each  metal  independently.  The  differentially  expressed  probe  sets 
were  selected  using  a  Benjamini-Hochberg  FDR  of  <0.001  and  a  fold  change  filter  >  1.8,  which 
identified  of  430,  456,  and  288  probe  sets  in  nickel,  chromium,  and  cadmium,  respectively 
(supplemental  data  table  S2  and  S3).  Many  of  the  probe  sets  are  differentially  expressed  in  more 
than  one  chemical  exposure,  yielding  a  total  of  992  differentially  expressed  probe  sets  (Figure  1). 

With  the  goal  of  comparing  the  response  to  these  metals  at  the  mechanistic  level,  a 
variety  of  enrichment  analyses  were  perfonned  to  identify  biological  processes  that  were 
statistically  over-represented  in  the  differentially  expressed  gene  lists.  Standard  enrichment 
analyses  were  perfonned  using  MetaCore  software  to  identify  transcription  factors  potentially 
associated  with  differentially  expressed  genes  in  our  data  set  (Figure  2)  and  Ingenuity  Pathway 
Analysis  software  to  identify  canonical  pathways  (Figure  3).  We  discovered  a  number  of 
enriched  transcription  factors  involved  in  DNA  damage  response,  cell  cycle,  cell  growth  and 
proliferation,  oxidative  stress,  and  hypoxia.  The  canonical  pathways  enriched  in  our  data  include 
processes  related  to  the  hypoxic  response,  glutathione  metabolism,  oxidative  stress,  and  retinoid 
receptor  signaling.  Upon  initial  review  of  these  findings  and  the  differentially  expressed  gene 
lists,  it  was  evident  that  several  biological  processes  are  represented,  but  conventional  gene 
ontology  and  pathway  categories  failed  to  capture  the  complexity  of  these  responses.  In  order  to 
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provide  a  more  comprehensive  view,  we  created  gene  ontology  bins  which  include  multiple  gene 
ontology  tenns  that  are  involved  in  the  same  biological  process  and  calculated  enrichment 
statistics  on  these  bins  (Table  1).  The  processes  that  are  enriched  in  our  data  set  include  oxidative 
stress,  DNA  damage,  apoptosis,  hypoxic  response,  and  energy  regulation.  These  are  consistent 
with  many  of  the  known  mechanisms  of  toxicity  for  the  three  metals  and  provide  a  starting  point 
to  allow  us  to  compare  and  contrast  the  response  among  the  metals. 

As  an  additional  step  in  categorizing  the  responses  to  the  toxicants,  we  perfonned  a 
cluster  analysis  across  all  three  of  the  metal  exposure  to  identify  potentially  co-regulated  genes 
using  Vxinsight  [15].  Three  clusters  were  identified  and  are  highlighted  in  white,  blue,  or  green 
(Figure  4),  and  contain  129,  456,  or  407  probe  sets  respectively.  The  blue  and  green  clusters  are 
primarily  comprised  of  probe  sets  that  are  up-  and  down-regulated,  respectively,  but  provide  no 
metal  specificity.  The  white  cluster  is  striking  since  it  is  tightly  clustered,  further  away  from  the 
others,  and  comprised  almost  entirely  of  probe  sets  up-regulated  in  response  to  chromium.  Many 
of  the  genes  in  this  cluster  are  involved  in  the  response  to  DNA  damage  as  discussed  below. 

Common  Response 

One  particular  interest  for  us  was  to  identify  processes  that  are  common  to  all  of  the 
metals.  In  these  experiments,  the  only  perturbed  biological  process  common  to  all  three  metals  is 
oxidative  stress,  a  known  effect  of  each  of  these  metals,  as  evidenced  by  the  oxidative  stress  bin 
being  enriched  in  response  to  the  three  metals  (Table  1)  [7].  While  some  of  the  changes  in  gene 
expression  are  consistent  across  the  metals,  our  observations  suggest  that  there  are  also  subtle 
variations  in  how  the  cells  respond  to  what  is  presumed  to  be  a  common  mechanism  of  toxicity. 
The  most  notable  differences  are  increased  production  of  the  anti-oxidant  protein  glutathione  in 
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response  to  nickel  and  ROS-induced  endoplasmic  reticulum  (ER)  stress  in  response  to  cadmium 
(Figure  5)  [16]. 

Oxidative  Stress 

The  nuclear  factor  erythroid  2 -related  factor  2  (Nrf2)  oxidative  stress  pathway  appears  to 
be  activated  in  response  to  all  three  metals.  The  pathway  is  statistically  enriched  only  for  nickel 
and  cadmium,  and  the  Nrf2  transcription  factor  is  statistically  identified  only  in  the  nickel  data. 
However,  key  Nrf2  controlled  genes  are  up-regulated  in  samples  from  all  three  metals,  including 
Hmoxl,  Sqstml,  and  glutathione-S-transferases.  Nrf2  is  a  transcription  factor  that  controls  the 
expression  of  important  detoxification  and  oxidative  stress  proteins  [17,18].  HMOX1  is  a 
ubiquitous  stress  response  protein  involved  in  reducing  the  effects  of  oxidative  stress  and 
apoptosis  [19].  Sequestosome  1  (SQSTM1)  has  been  shown  to  play  a  role  in  the  sustained 
activation  of  Nrf2  in  response  to  oxidative  stress  [20]. 

Other  genes  and  transcription  factors  known  to  respond  to  oxidative  stress  outside  of  the 
Nrf2  pathway  are  also  differentially  regulated  by  all  three  metals.  For  example,  Haol  is  down- 
regulated.  The  enriched  transcription  factors  include  FKHR  and  NF1C  for  cadmium,  HSF1  and 
ATF-4  for  both  nickel  and  cadmium,  and  OCT1  for  all  three  chemicals.  HAOl  is  a  liver  specific 
enzyme  that  converts  a-hydroxy  acids  to  a-keto  acids  while  reducing  molecular  oxygen  to  H2O2, 
and  has  been  shown  to  be  down-regulated  due  to  oxidative  stress  [21].  FKHR  has  been  shown  to 
be  a  principal  component  in  the  response  to  oxidative  stress  by  stimulating  the  expression  of 
metal  containing  antioxidant  proteins  [22].  HSF1  decreases  intracellular  reactive  oxygen  species 
generation,  thereby  protecting  against  further  damage  [23]. 

The  enrichment  of  the  Nrf2-mediated  oxidative  stress  response  canonical  pathway  and 
the  modulation  of  key  genes  known  to  respond  to  oxidative  stress  suggest  that  all  three  metals 
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induce  oxidative  stress,  with  chromium  showing  the  lowest  level  of  induction.  The  metals  do, 
however,  differentially  affect  other  mechanism  that  control  oxidative  stress 

Glutathione  pathways  appear  to  be  activated  only  in  response  to  nickel.  The  IPA 
Glutathione  Metabolism  canonical  pathway  (Figure  3)  is  significantly  enriched  due  only  to 
exposure  to  nickel,  and  a  number  of  the  genes  affected  solely  by  nickel  are  involved  in  the 
protection  of  the  cell  by  glutathione  (Fig  4),  including  glutamate-cysteine  ligase,  modifier 
subunit  (Gclm)  and  ATP-binding  cassette  sub-family  C  member  1  {Abed).  GCLM  is  the  first 
and  rate  limiting  enzyme  of  glutathione  synthesis  and  ABCC 1  has  been  shown  to  be  a 
glutathione  transporter  [24,25].  The  role  of  glutathione  in  response  to  nickel  toxicity  is  likely 
two-fold;  both  as  an  antioxidant  and  in  neutralizing  the  toxic  effects  of  nickel  by  acting  as  a 
chelator,  thus  increasing  the  efflux  of  nickel  out  of  the  cell  [26].  This  up-regulation  of  glutathione 
metabolism  suggests  a  mechanism  unique  to  nickel. 

Oxidative  stress  in  cadmium  exposed  cells  leads  to  ER  stress  including  the  induction  of 
the  unfolded  protein  response  and  apoptosis.  A  number  of  genes  modulated  solely  by  cadmium 
are  indicative  of  the  unfolded  protein  response,  a  mechanism  not  seen  in  response  to  nickel  or 
chromium  (Figure  5).  Six  genes  encoding  chaperones  ( Hspala ,  Hspalb ,  Hspb8,  Dnajbl, 

Dnajc3,  and  Cryab )  are  up-regulated.  Chaperone  proteins  are  known  to  be  involved  in  apoptosis 
as  well  as  the  folding  and  degradation  of  damaged  proteins  in  the  unfolded  protein  response  [27]. 

ER  stress  can  lead  to  apoptosis,  and  our  data  support  the  presence  of  apoptosis  in 
cadmium  exposed  cells.  The  apoptosis  ontology  bin  is  significantly  enriched  (Table  1),  and 
Caspase  4  ( Casp4 ),  encoding  an  apoptosis-related  cysteine  peptidase  [28],  is  up-regulated  due  to 
exposure  to  cadmium  (Figure  5).  Two  other  apoptotic  genes,  typically  seen  up-regulated  in 
response  to  DNA  damage,  were  also  up-regulated  in  the  cadmium  data:  a  protein  phosphatase  1 
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regulatory  (inhibitor)  subunit  15A  ( Ppplrl5a ),  and  DNA-damage  inducible  transcript  3  ( Ddit3 ) 
[29,30].  However,  in  this  work  we  believe  that  the  induction  of  Ppplrl5a  and  Ddit3  is  purely 
related  to  their  role  in  apoptosis  and  not  indicative  of  DNA  damage.  Overall,  these  results 
suggest  that  cadmium-induced  oxidative  stress  causes  ER  stress  leading  to  the  unfolded  protein 
response  and  apoptosis. 

While  the  three  metals  do  share  some  common  responses  to  the  disturbance  of  the  cell’s 
nonnal  redox  state,  each  metal  affects  a  unique  subset  of  genes.  Chromium  appears  to  have  a 
lower  level  of  induction  for  the  Nrf2  pathway,  nickel  induces  production  of  the  antioxidant 
glutathione,  and  cadmium  mediates  an  RO-  induced  ER  stress  characterized  by  the  unfolded 
protein  response  and  apoptosis. 

Unique  Responses 

In  addition  to  shared  responses,  we  were  interested  in  identifying  mechanisms  of  toxicity 
that  are  unique  to  each  metal.  Based  on  the  gene  expression  changes  present  in  our  data, 
chromium  is  unique  in  causing  DNA  damage;  nickel  causes  a  hypoxic  response  and  perhaps 
disruption  of  chromatin  structure;  and  cadmium  causes  a  disruption  of  retinoic  acid  signaling 
pathways. 

Cr-induced  DNA  Damage 

Chromium  is  the  only  metal  of  the  three  that  appears  to  cause  DNA  damage.  Our  data 
supports  this  mechanism  with  genes  involved  in  DNA  repair  and  replication  being  modulated  by 
exposure  to  the  chromium  (Figure  6)  and  the  DNA  damage  ontology  bin  being  enriched  (Table 
1).  The  transcription  factor  E2F1,  which  is  induced  by  DNA  damage,  plays  an  important  role  in 
DNA  repair  at  stalled  replication  forks  [31].  Strikingly,  among  the  probe  sets  in  the  white 
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Vxlnsight  cluster  (Figure  4),  most  of  which  are  induced  by  chromium,  40  of  the  79  probe  sets 
with  annotation  provided  by  Affymetrix  are  involved  in  the  DNA  damage  response.  Many  of 
these  genes  have  products  that  can  be  found  in  the  DNA  synthesome,  which  is  a  multiprotein 
complex  involved  in  DNA  replication  [32,33].  The  genes  proliferation  cell  nuclear  antigen 
(Pena),  replication  protein  A  ( Rpa ),  the  minichromosome  maintenance  complex  component 
genes,  which  encode  helicases,  DNA  ligase  (Ligl),UNK  polymerase  s  (Pole),  and  DNA 
polymerases  6  ( Pole l)  are  all  up-regulated  as  a  result  of  exposure  to  chromium,  suggesting  that 
there  is  an  increase  in  DNA  synthesis,  likely  due  to  chromium-induced  DNA  damage. 

These  proteins  that  make  up  part  of  the  synthesome  play  specific  roles  in  DNA  synthesis 
and  repair  that  are  consistent  with  the  known  mechanism  of  chromium-induced  DNA  damage. 
Chromium  has  been  shown  to  directly  interact  with  DNA  and  cause  damage,  by  fonning  DNA 
adducts  and  causing  DNA  strand  breaks  [7,34].  DNA  polymerases  5  and  s  are  associated  with 
proof-reading  and  repair  activity  [35].  These  polymerases,  as  well  as  several  other  proteins 
including  RPA2  and  PCNA,  may  be  involved  in  excision  repair  to  remove  DNA  adducts.  DNA 
ligase  is  involved  in  repairing  double  strand  breaks,  which  are  known  to  accumulate  due  to 
chromium  toxicity.  Some  of  the  encoded  proteins  of  the  genes  differentially  expressed  in  the 
DNA  synthesome,  such  as  those  fonning  the  minichromosome  maintenance  complex  (MCM) 
and  origin  recognition  complex  (ORC),  are  involved  specifically  in  the  initiation  of  DNA 
synthesis  [36].  Since  it  has  been  shown  that  chromium-DNA  adducts  reduce  the  ability  of 
synthesome  to  initiate  replication  [37],  the  expression  of  these  genes  may  be  up-regulated  as  the 
cell  attempts  to  repair  damage  caused  by  the  chromium.  These  responses  suggest  that  the  direct 
interaction  of  chromium  with  DNA  and  the  formation  of  adducts  and  stand  breaks  are  involved 
in  chromium  toxicity. 
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Although  nickel,  cadmium,  and  chromium  believed  to  be  geno toxic,  we  did  not  detect 
changes  in  gene  expression  attributable  to  DNA-damage  in  cadmium  and  nickel  exposed  cells 
despite  having  clear  evidence  for  ongoing  DNA  replication  and  repair  caused  by  chromium.  At 
the  concentrations  tested,  chromium  is  the  only  of  the  three  metals  to  cause  high  levels  of  DNA 
damage  in  H-II-4-E-C  cells. 

Ni/Fe  substitution 

Gene  expression  changes  seen  in  the  nickel-exposed  samples  were  consistent  with  a 
hypoxic  response.  The  hypoxic  response  ontology  bin  (Table  1),  HIF-la  canonical  pathway 
(Figure  3),  and  HIF-la  transcription  factor  (Figure  2)  are  enriched  in  the  nickel  data.  HIF-la  is  a 
transcription  factor  which  induces  the  transcription  of  genes  involved  in  glycolysis,  glucose 
transport,  apoptosis,  and  other  cellular  process  as  a  result  of  a  change  in  the  intracellular  oxygen 
concentration  (Beyersmann  and  Hartwig  2008).  Additionally,  the  glycolysis/gluconeogeneis 
canonical  pathway  and  energy  regulation  ontology  bins  were  both  enriched  (Table  1).  These 
could  also  be  potentially  a  result  of  HIF-la  regulation,  as  hypoxic  conditions  and  HIF-la 
activation  are  known  to  interfere  with  cellular  energy  metabolism  such  as  glycolysis,  causing  a 
cell  to  shift  toward  nonoxidative  forms  of  ATP  production  and  enhancing  production  of 
glycolytic  enzymes  and  glucose  transporters  [38].  Furthennore,  the  genes  encoding  lactate 
dehydrogenase  A  ( Lclha ),  pyruvate  dehydrogenase  ( Pdkl ),  phosphoglycerate  kinase  1  ( Pgkl )  and 
solute  carrier  family  2  (facilitated  glucose  transporter)  member  1  ( Slc2al )  are  up-regulated;  all 
which  contribute  to  low  oxygen  energy  utilization  (Figure  5)  [38].  These  data  suggests  that  nickel 
alters  the  expression  of  known  HIF-la  targets  and  induces  a  hypoxic  response.  The  cause  of  the 
hypoxic  response  in  the  case  of  nickel  exposure  may  not  be  low  oxygen  levels.  It  is  thought  that 
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nickel  activates  HIF-la  by  preventing  the  degradation  of  the  protein  either  through  the  depletion 
of  ascorbate  or  by  replacing  iron  in  the  hydroxylases  responsible  for  HIF-la  degradation  [39,40]. 

A  number  of  the  genes  up-regulated  specifically  by  nickel  are  involved  in  chromatin 
structure  modifications,  including  two  jumonji  family  histone  demethylases  ( Jmjclla  and  Jmjd6). 
It  has  been  shown  that  nickel  may  inactivate  jumonji  family  histone  demethylases  by  replacing 
iron  in  the  enzyme’s  active  site,  and  the  increase  seen  may  be  due  to  this  inactivation  and  not 
DNA  damage  [41].  This  is  a  second  example  of  nickel  substituting  for  iron  as  a  cofactor  that 
disrupts  nonnal  function  of  a  protein. 

Chromatin  structure  and  the  hypoxic  response  are  affected  by  nickel  exposure  alone. 
While  these  two  functions  are  clearly  unrelated,  they  are  both  mediated  by  enzymes  that  require 
iron  as  a  cofactor.  It  is  known  that  nickel  can  substitute  for  iron  in  many  enzymes  and  block  their 
function  [42].  In  nickel  exposed  cells,  Ni/Fe  substitution  may  instigate  the  induction  of  the 
hypoxic  response  and  changes  in  expression  of  genes  related  to  chromatin  structure.  Nickel’s 
ability  to  replace  iron  as  a  cofactor  in  multiple  proteins  is  a  major  contributor  to  nickel  toxicity 
that  is  unique  amongst  the  metals  in  this  study. 


Retinoic  acid  signaling 

The  data  also  suggests  a  mechanism  unique  to  cadmium:  retinoic  acid  signaling.  The  IPA 
canonical  pathways  FXR/RXR  Activation  and  RXR  Activation  were  significant  only  in  response 
to  cadmium  (Figure  3).  Transcription  factors  comprised  of  the  retinoid  family  receptors 
including  RAR  gamma,  RAR  alpha,  ROR  alpha,  and  the  vitamin  D  receptor  (VDR)  are  enriched 
in  response  to  only  cadmium  (Figure  2).  Retinoic  acid  is  a  hormone-like  molecule  that  is 
involved  in  the  regulation  of  cell  differentiation  and  proliferation  whose  effects  are  mediated  by 
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retinoic  acid  receptors  [43],  It  has  been  suggested  that  cadmium  acts  as  an  environmental 
teratogen  by  increasing  the  amount  of  retinoic  acid  through  interference  with  the  retinoic  acid 
metabolizing  genes  [44].  These  enriched  pathways  and  transcription  factors  suggest  that  the 
disruption  of  retinoid  family  signaling  is  a  cadmium  specific  mechanism. 

Unexpected  Findings 

While  we  were  able  to  identify  both  common  and  unique  responses  to  nickel,  cadmium, 
and  chromium,  each  metal  also  has  known  mechanisms  that  we  expected  our  data  to  reflect 
based  on  the  literature,  but  were  not  apparent.  We  did  not  find  evidence  of  nickel-  and  cadmium- 
induced  DNA  damage  or  a  strong  induction  of  oxidative  stress  by  chromium,  all  of  which  are 
well  documented  effects  in  other  systems. 

Nickel,  cadmium,  and  chromium  have  all  been  shown  to  cause  to  DNA  damage,  in  our 
data,  however,  only  chromium  was  genotoxic.  The  DNA  ontology  bin  (Table  1)  and  E2F1 
transcription  factor  (Figure  2)  are  significant  only  for  chromium,  and  genes  involved  in  DNA 
metabolism  were  up-regulated  only  in  response  to  chromium  (Figure  6).  Chromium  can  directly 
interact  with  and  damage  DNA,  while  nickel  and  cadmium  only  indirectly  damage  DNA  through 
the  fonnation  of  reactive  oxygen  species  and  by  interfering  with  DNA  repair  enzymes.  Our 
inability  to  detect  the  nickel  and  cadmium  induced  DNA  damage  may  be  a  result  of  these 
mechanisms;  a  24  hour  exposure  period  may  not  have  been  long  enough  for  DNA  damage  to 
accumulate  in  H— II-E-C3. 

We  also  did  not  observe  evidence  of  a  strong  induction  of  oxidative  stress  due  to 
exposure  to  chromium.  While  the  Nrf2  oxidative  stress  response  is  enriched  in  the  cadmium  and 
nickel  data  (Table  1;  Figure  5),  and  the  NrfZ  transcription  factor  is  enriched  in  response  to  nickel 
(Figure  2),  neither  are  significantly  enriched  in  response  to  chromium.  Also,  the  change  in 
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magnitude  of  some  of  the  key  genes  involved  in  the  response  to  oxidative  stress  is  not  as  large 
due  to  chromium  as  it  is  for  nickel  and/or  cadmium.  For  example,  Hmox  is  increased  almost  7 
and  4.5  fold  in  response  to  nickel  and  cadmium  (high  dose),  respectively,  but  only  2  fold  due  to 
chromium  (high  dose)  (Figure  5).  The  lack  of  evidence  supporting  oxidative  stress  due  to 
chromium  as  compared  to  nickel  and  cadmium  is  surprising  as  the  fonnation  of  reactive  oxygen 
species  as  Cr(VI)  is  reduced  to  Cr(III)  intracellularly  is  well  documented  [45].  Perhaps  at  the 
concentrations  used  in  this  work,  chromium-induced  DNA  damage  was  the  dominant  effect  of 
the  metal,  thus  overshadowing  the  oxidative  stress  response,  or  the  response  may  be  unique  to 
H4-II-E-C3  cells. 

Caveats  of  analysis 

In  evaluating  the  results  of  our  analysis,  there  are  several  important  worth  noting.  The 
first  is  that  the  number  of  enriched  categories  appears  somewhat  smaller  in  the  chromium  data 
set,  which  we  believe  might  be  due  in  part  to  a  weakness  in  the  enrichment  analysis  approach. 
The  second  is  that  the  high  overlap  in  genes  across  many  pathways  may  lead  to  the  statistical 
enrichment  of  processes  which  are  truly  uninvolved.  The  final  is  the  need  for  equipotent 
concentrations  across  the  metals  to  allow  a  realistic  comparison  of  the  toxic  mechanisms. 

In  the  IPA  canonical  pathway  analysis  for  chromium,  only  4  pathways  were  significant 
compared  to  9  and  24  from  cadmium  and  nickel,  respectively,  and  we  observed  a  paucity  of 
enriched  transcription  factors  due  to  chromium  exposure  as  compared  to  the  nickel  and  cadmium 
exposure.  Since  a  large  proportion  of  the  chromium  modulated  genes  are  involved  in  the  DNA 
damage  response,  it  may  have  dominated  the  enrichment  analysis,  masking  other  biologically 
important  perturbed  processes.  If  the  DNA  damage  genes  are  removed  from  the  chromium 
analysis,  additional  canonical  pathways  achieve  statistical  significance,  including  Glycerolipid 
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Metabolism,  Glycolysis,  Starch  and  Sucrose  Metabolism,  and  FXR/RXR  Metabolism  (data  not 
shown). 

Moreover,  applying  strict  criteria  for  differential  expression  as  was  done  in  this  work  can 
restrict  the  number  of  enriched  pathways.  Less  stringent  criteria  for  differential  expression  might 
have  increased  the  number  of  enriched  pathways  by  increasing  the  number  of  genes  contributing 
to  the  enrichment  analysis.  The  stricter  criteria  for  differential  expression  used  here  may  have  led 
to  a  high  false  negative  rate,  but  our  findings  are  well  supported. 

The  large  overlap  existing  among  the  genes  in  different  pathways  and  among  the 
transcription  factor  target  lists  can  also  complicate  the  interpretation  of  enrichment  analyses.  The 
jun  proto-oncogene  {Jun),  considered  a  “hub”  molecule,  is  present  in  85  IPA  canonical 
pathways,  and  Hmox  is  present  in  12.  In  the  transcription  factor  enrichment  analysis  factors  with 
similar  binding  sites,  such  as  USF1  and  USF2,  and  the  retinoic  acid  receptors  RAR-gamma, 
RAR-alpha,  and  ROR-alpha,  are  all  enriched.  It  is  possible  that  several  of  these  transcription 
factors  may  have  been  assigned  to  an  individual  gene  because  of  a  single  binding  site.  A  few 
differentially  expressed  genes  could  therefore  cause  significant  enrichment  of  many  different 
pathways  or  transcription  factors,  leading  to  an  incorrect  analysis.  Thus,  infonnation  gained  from 
simple  enrichment  analysis  tools  must  be  viewed  with  caution.  To  prevent  the  inclusion  of 
spurious  processes,  we  have  analyzed  individual  genes  within  the  enriched  categories  to  help 
ensure  that  the  results  reported  here  are  biologically  relevant. 

An  extremely  important,  yet  challenging,  component  of  performing  a  comparative 
toxicogenomic  analysis  is  setting  equipotent  stimuli  across  the  study  conditions.  In  the  toxicity 
study  reported  here,  we  found  that  the  viability  assays  used  in  our  range  finding  did  not  show 
good  concordance  with,  or  similar  sensitivity  to,  the  measurement  of  gene  transcripts. 
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Additionally,  each  metal  affected  the  cells  differently,  and  concentrations  of  the  metals  at  equal 
levels  of  cytotoxicity  did  have  the  same  effect  at  the  gene  expression  or  biological  process  level. 
The  same  biological  processes  could  be  perturbed  at  different  levels  of  cytotoxicity  for  the 
different  metals.  At  the  concentrations  we  tested,  similar  numbers  of  genes  were  differentially 
expressed  among  the  metal  exposures.  Therefore,  based  on  the  similar  number  of  differentially 
expressed  genes  and  the  similar  levels  of  cytotoxicity,  we  believe  that  we  approximated 
equipotency  sufficiently  well  to  produce  useful  results. 

Conclusion 

Nickel,  chromium,  and  cadmium  are  heavy  metals  commonly  found  in  industry  use  and 
in  the  environment  that  have  adverse  health  effects.  In  order  to  identify  common  and  unique 
molecular  mechanisms  of  toxicity  for  each  metal,  a  microarray  study  was  performed  using  rat 
hepatoma-derived  cells  exposed  to  the  metals.  Nickel,  cadmium,  and  chromium  all  induced 
common  effects  when  broadly  viewed,  but  the  detailed  mechanisms  and  pathways  involved  were 
unique  to  the  metals.  All  three  metals  cause  oxidative  stress,  and  the  cells  response  to  it  was 
mediated  at  least  in  part  through  the  Nr  IT  transcription  factor.  However,  the  oxidative  stress 
response  was  distinct  for  each  metal.  Chromium  had  the  lowest  level  of  response,  nickel  induced 
synthesis  of  the  anti-oxidant  glutathione,  and  cadmium  lead  to  ROS  mediated  ER  stress  and  the 
unfolded  protein  response.  Further,  all  three  metals  are  known  to  be  genotoxic,  yet  in  this  work, 
only  chromium  caused  extensively  stimulation  of  DNA  repair  mechanisms,  likely  through  DNA 
adduct  fonnation  and  DNA  strand  breakage.  Substitution  of  nickel  in  place  of  iron  as  an  enzyme 
co-factor  caused  Hif-la  activation  and  disruption  of  chromatin  structural  proteins,  which  was  a 
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mechanism  unique  to  metal.  Cadmium  caused  disruption  of  retinoic  acid  signaling,  which  is  a 
likely  mechanism  for  cadmium-induced  teratogenicity. 

In  conclusion,  the  gene  expression  of  the  H4-II-E-C3  cell  line  was  investigated  to  identify 
mechanism  of  toxicity  for  nickel,  chromium,  and  cadmium.  Identified  mechanisms  included 
oxidative  stress,  DNA  damage,  Ni/Fe  substitution,  and  disruption  of  retinoic  acid  signaling. 

While  microarray  analysis  suggests  toxicity  pathways  involved  in  heavy  metal  intoxication, 
further  investigation  will  be  required  to  verify  these  findings.  This  work  provides  a  starting  point 
for  future  studies  by  provide  key  genes  and  transcription  factors  that  may  be  directing  the  cells’ 
response  to  toxic  insults  by  nickel,  chromium,  or  cadmium. 
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